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The finite temperature transition or crossover in QCD witii two degenerate Wilson quarks is investigated using 
a renormalization group improved action. At /? = 2.0 and 2.1 where ~ 1.0 — 1.2 GeV, the expectation value 
of the Polyakov loop and the pion screening mass on an 8'^ x 4 lattice vary smoothly with the hopping parameter 
through the transition/crossover. The quark screening mass in the high temperature phase agrees well with that 
in the low temperature phase calculated on an 8'* lattice. The smooth transition of the observables is totally 
different from the sharp transition found for the standard action at /3 = 5.0 and 5.1 where is also 1.0 — 1.2 
GeV. 



1. Introduction 

Recently MILC group has extensively investi- 
gated the finite temperature transition/crossover 
in the case of two degenerate Wilson fermions and 
found unexpected phenomena0J^. For Nt—A, 
when the quark mass is heavy the transition is 
smooth and when the quark mass is light it is 
also smooth. However, in the intermediate mass 
region(/3 5.0), the transition is very sharp. This 
is completely opposite to what was supposed to 
be realized: It was supposed to become weaker 
or disappear. Furthermore, for Nt — 6, in the 
range of intermediate mass, they observed clear 
two state signals. 

We reproduce in fig.^-a the data[0 for the pion 
screening mass squared m^, the expectation value 
of the Polyakov loop L and the quark screening 
mass niq at /3 = 5.0, together with our data|^. 
The transition of the Polyakov loop is sharp and 
a cusp appears in the pion screening mass at the 
transition point. The quark mass defined through 
an axial- vector Ward identity |^,^ also exhibits 
peculiar behavior. The ruq in the high temper- 
ature phase on the Nt — 4 lattice behaves singu- 
larly and does not agree with that in the low tem- 
perature phase. This is in clear contrast with the 
fact that the rUq is almost independent of whether 
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the system is in the high temperature phase or in 
the low temperature phase at /3 = 5.85 in the 
quenched QCD § and at /3 = 5.5 for Nf = 2 
caseQ. 

One possible reason for these phenomena is 
that we are far from the continuum limit. In 
this work, we study the finite temperature tran- 
sition/crossover in the case of two degenerate 
quarks using a renormalization group (RG) im- 
proved action where we are supposed to be close 
to the continuum limit. Our main concern is to 
see whether the unexpected phenomena described 
above persist with the improved action. 

2. Improved action 

The RG improved action for pure gauge theory 
we proposed about ten years ago is 

Sg = l/.g^{co y^(plaquette) + ci ^(1 x 2 loop)} 

with ci — —0.331 and cq — 1 — 8ci. Each oriented 
loop appears once in the sum. The form of the 
gauge action is determined by a block-spin RG 
analysis in the pure gauge theory [||. We calcu- 
lated the string tension [p|jic|] and the hadron spec- 
trumflprl and investigated the topological prop- 
erties and the C/(l) problem]!^ in the quenched 
QCD with this improved action. We observe the 
following among others: The approach to the per- 
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Figure 1. Observables on iVt = 4 (filled symbols) and iVt = 6 or 8 (open symbols) lattice, a) At /3 = 5.0 
with the standard action. Data are taken from 0| and b) At /3 = 2.0 with the improved action. 



turbative result for the expectation value of the 
plaquette is much faster than the case of the stan- 
dard one-plaquette action Also note that the 
scale parameter of the improved action Kim is 
close to the k-^: Ajj^/Ajm ^ 0.488@. There- 
fore the bare coupling of the improved action is 
already an improved coupling jl^. 

With these results it may be worthwhile to use 
a RG improved action for full QCD. The action 
we use as an improved action for full QCD is given 

by 

o nII\I I cW't^son 

'^gauge ' ^quark ' 

The action is a sum of the RG improved action 
for gluons and the Wilson action for quarks. This 
action reduces to the improved pure gauge action 
in the heavy quark mass limit. There is also a 
possibility to improve the quark action. However 
we take the Wilson action for the quark action 
as a first step, because we believe that the ef- 



fect of the improvement of the gauge sector is 
much more significant than that of the quark sec- 
tor by the following reasoning: After integration 
over fermion variables, the quark action may be 
written as a sum over various types of Wilson 
loops, assuming the convergence of hopping pa- 
rameter expansion. The coefficients of various 
types of Wilson loops thus obtained are much 
smaller compared with ci = —0.331 for the im- 
proved pure gauge action. 

3. Numerical Simulations 

Numerical simulations are done at /3 = 2.0 and 
2.1 on an 8'^ X 4 lattice at several hopping pa- 
rameters. In order to fix the scale and to com- 
pare the results in the high temperature phase 
with those in the low temperature phase, we also 
make simulations on an 8^ x 8 lattice. We use the 
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hybrid Monte Carlo algorithm to generate gauge 
configurations with molecular dynamics step size 
Ar — 0.01. Momentum refresh is done at every 
one unit of simulation time. Acceptance rates are 
larger than 0.92 for all cases. We discard first 100 
^ 400 trajectories for thermalization and use 100 
~ 1000 trajectories for measurement. The num- 
bers of trajectories for thermalization and mea- 
surement are chosen taking into account whether 
the simulation point is close or not to the transi- 
tion point. Wilson loops and Polyakov loops are 
measured at every trajectory. Screening correla- 
tion functions are measured at every 5 trajectories 
on an 8^ x 16 x Ni lattice obtained by doubling 
the original lattice. 

Using the p meson mass calculated on the 
S'^ X 16 lattice as an input, we obtain a"^ = 1.01 
GeV and 1.18 GeV at /? = 2.0 and 2.1, respec- 
tively. Therefore these /?'s correspond approxi- 
mately to /3 — 5.0 and 5.1 of the standard action, 
respectively. 

4. Results 

Fig.^b shows the results for the observables 
at /? = 2.0 versus XjK — 1/Kc, where the Kc is 
determined from the quark mass on the Nt = A 
lattice. The data indicate the following: 1) The 
transition/crossover of the Polyakov loop around 
l/K - l/Kc ~ 0.5 is very smooth. 2) The hop- 
ping parameter dependence of the quark mass is 
also smooth and its value on the A^t = 4 lattice 
agrees well with that on the Nt — 8 lattice. This 
implies that the value of the quark mass is al- 
most independent of whether the system is in the 
low temperature phase or in the high tempera- 
ture phase. 3) The pion screening mass does not 
exhibit a cusp. 

These features of the transition/crossover at 
(3 = 2.0 are totally different from those observed 
for the standard action at (3 — 5.0. Note that 
are about 1 GeV for both cases. We also observe 
a similar very clear contrast between the results 
at /? = 2.1 with the improved action and those at 
/3 = 5.1 with the standard action where are 
both approximately 1.2 GeV. 

Furthermore preliminary results in the range of 
of/? = 1.6 — 2.15 indicate that the transitions are 



smooth from /3 = 2.15 down to at least /3 = 1.8. 
Therefore we conclude that the sharp transition 
in the intermediate mass region disappears with 
the improved action at Nt = 4. This suggests a 
strong possibility that the first order phase transi- 
tion observed by MILC at TVj = 6 disappears with 
this action. To confirm this, we are planning to 
extend the work to a wider range of parameter, 
in particular, on a lattice with Nt = 6. 

The simulations are performed with Fujitsu 
VPP500/30 at University of Tsukuba. We would 
like to thank C. DeTar for providing us with data 
by MILC collaboration. This work is in part sup- 
ported by the Grant-in- Aid of Ministry of Edu- 
cation, Science and Culture (N0.O6NPO6OI). 
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